home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / DDJMAG / DDJ9203.ZIP / 80X87.ZIP / TESTQUAD.C < prev    next >
Text File  |  1991-11-09  |  878b  |  31 lines
  1. /* testquad.c: Program to test the correctness of the quadratic-
  2.  *           solving routine in 80x87 assembler, and to illustrate
  3.  *           its use.
  4.  * Copyright (C) 1991 by Nicholas Wilt.  All rights reserved.
  5.  */
  6.  
  7. #include "testfpu.h"
  8.  
  9. void
  10. main()
  11. {
  12.   char inp[256];    /* Array for input string */
  13.  
  14.   int done = 0;     /* Termination var */
  15.   double a, b, c;    /* Coefficients of quadratic */
  16.   double x1, x2;    /* Roots of quadratic, if any */
  17.  
  18.   while (! done) {
  19.     printf("a = ");    gets(inp);    a = atof(inp);
  20.     printf("b = ");    gets(inp);    b = atof(inp);
  21.     printf("c = ");    gets(inp);    c = atof(inp);
  22.     if (solve_quadratic(a, b, c, &x1, &x2))
  23.       printf("Roots are %.2f and %.2f\n", x1, x2);
  24.     else
  25.       printf("No real roots.\n");
  26.     printf("Test another (Y/N)? ");
  27.     gets(inp);
  28.     done = ! (toupper(inp[0]) == 'Y');
  29.   }
  30. }
  31.